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Aims. The availability of wide angle and deep surveys, both in the optical and the radio band, allows us to explore the evolution of radio sources 
with optical counterparts up to redshift z ~ 1.1 in an unbiased way, using not only large numbers of radio sources but also well defined control 
samples of radio quiet objects. 

> Methods. We use the 1.4 GHz VIMOS-VLA Deep Survey and the optical VIMOS-VLT Deep Survey and the CFHT Legacy Survey to compare 
^ ' the properties of radio loud galaxies with respect to the whole population of optical galaxies. The availability of multiband photometry and high 
I I ' quality photometric redshifts allows to derive rest frame colors and radio luminosity functions down to a limit of a B rest-frame magnitude of 
, Mb = -20. Moreover, we derive spectrophotometric types, following the classification of Zucca et al. (2006), in order to have a priori knowledge 
T-H • of the optical evolution of different galaxy classes. 
^ Results. Galaxy properties and luminosity functions are estimated up to z ~ 1 for radio loud and radio quiet early and late type galaxies. Radio 
loud late type galaxies show significantly redder colors with respect to radio quiet objects of the same class and this is an effect related to the 
presence of more dust in stronger star forming galaxies. Moreover, we estimate optical luminosity functions, stellar masses and star formation rate 
distributions for radio sources and compare them with those derived for a well defined control sample, finding that the probability for a galaxy 
T-H to be a radio emitter significantly increases at high values of these parameters. Radio loud early type galaxies show luminosity evolution of their 
• bivariate radio-optical luminosity function, due to an evolution in the radio-optical ratio. The lack of evolution of the mass function of radio loud 
early type galaxies means that no new AGN are formed at redshift z < 1. On the contrary, radio loud late type objects show a strong evolution, 
both in luminosity and in density, of the radio luminosity function for z > 0.7. This evolution is the direct effect of the strong optical evolution of 
this class and no significant change with redshift of the radio-optical ratio is required. With the knowledge of the radio-optical ratio and the optical 
and radio luminosity functions for late type galaxies, we show that it is possible to estimate the star formation history of the Universe up to redshift 
, z ~ 1.5, using optical galaxies as tracers of the global radio emission. 

Key words. Galaxies:fundamental parameters - Galaxies:general - Galaxies:luminosity function - Radio continuum: galaxies 

1 . Introduction providing different results (ISmolcic et al.L l2008l and references 

therein). 

Deep 1.4 GHz counts show an upturn below a few milli- 

jansky (mJy), coiTesp onding to a rap id increase in the num- I" f^'^^' *e photometric and spectroscopic work needed in 

ber of faint sources (IWindhorst et all [T985h . First results of *e optical identification process is very demanding in terms of 

the spectroscopic follow-up of relatively bright counterparts telescope time, since a significant fraction of faint radio sources 

(B < 22) identified many of them as blue g alaxies with spectra ^as also very faint optical counterparts. It is therefore clear that 

showi ng evidence of intense star formation (F ranceschini et aH o^'^er to investigate the nature and evolution of the sub-mJy 

[19881 iBenn etaLL [T993h . But with increasing depth of the Population it is necessary to couple deep radio and optical (both 

optical follow-up, it resulted that an increasing number of imaging and spectroscopic) observations over a reasonably large 

sub-mJy sources was associat ed with earlie r type galaxies °^ the sky. 

(Gruppioni et al., 2003; Prandoiiietin 120011: iHammer et"an The properties of the radio s ources in the loca l Universe 

ri995HAfonso et al.L 120061) . have been extensively studied bv iBest et al] (|2005|) using the 

The relative fractions of the various populations responsi- NVSS/FIRST radio surv eys coupled with the opti cal Sloan 

ble of the sub-mJy radio counts (AGN, starburst, late and early Digital Sky Survey and bv lMagliocchetti et all (|2002|) using the 

type galaxies) are still amatter of debate, with diff^erent surveys same radio data and the 2dF galaxy redshift survey for the op- 
tical part. These optical surveys permit a spectral classification 

for the optical counterparts of radio sources, mainly dividing the 

Send offprint requests to: S.Bardelli sandro.bardelli@oabo.inaf it sample in early and late type galaxies (broadly representative of 
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passive and star-forming galaxies) on the basis of spectral fea- 
tures such as the 4000 A break or colors. The results, represented 
by the bivariate radio-optical luminosity functions, show that the 
luminosity function of early type galaxies (presumably hosting 
an AGN) can be parametrized with a power law function flatter 
than that of late type galaxies. 

Moving to higher redshift, the Phoenix Deep survey al- 
lowed the study of th e evolution of the radio luminosity function 
jAfonso et al.L l2005l) and showed that the radio emission of star- 
forming galaxies strongly evolved from z ~ 0.1 up to z ~ 0.5. 
More recent surveys explored th e star formation density evolu- 
tion at even higher re dshifts (Smolcic et al., 2008; Tvison et al.L 
120071: ISevmour et all [2008; Barger et al.., .2007.) finding the typ- 
ical trend of rapidly increasing star formation from redshift zero 
to ~ 0.8 - 1 .0 and then a plateau (or a decrease) for higher red- 
shifts. 

When one wants to study radio sources as a function of red- 
shift, it is necessary to go through the radio-optical identifica- 
tion process. For this reason, it is critical to know as accurately 
as possible the properties of the optical counterparts in order to 
discriminate which is an effect due to the optical or to radio evo- 
lution. 

The VIMOS-VLT Deep Survey (iLe Fevre et al.Ll2005h . com- 
plemented with photometric data from the CFHT Legacy 
Survey, gives a unique opportunity to study the evolution of opti- 
cal galaxies. In fact, it offers the availability of a large sample of 
reference redshifts obtained from a deep spectroscopic survey, 
a large set of photometric bands and high quality photometric 
redshifts over a wide area. 

This allows us to study the properties of the coun t erparts 
of an assoc i ated 1 .4 GHz VLA survey (iBondi et al.L 120031: 
ICiheei et al.L llool using the rest-frame colors and therefore 
avoiding the use of redshift dependent quantities. 

The aim of this paper is to study the optical properties of 
galaxies hosting a radio source with respect to a well defined 
control sample of galaxies, for which the selection effects and 
the evolutionary behaviour are known a priori, in order to un- 
derstand the relations between the radio and optical properties. 
For this reason, in this paper we define as "radio galaxies" or 
"radio loud galaxies" all galaxies with detected radio emission: 
therefore our definition is based only on observed quantities. 

The advantage of this work is the simultaneous availability 
of a deep radio survey, a multiwavelength optical band coverage 
over a relatively large area, accurate photometric redshifts and 
absolute magnitudes which minimize the model dependency, but 
the most important point is the availability of a control sample 
with well studied luminosity and stellar mass functions as a func- 
tion of redshift and galaxy type, from which it is possible to ex- 
tract in a completely homogeneus way the galaxies which host a 
radio source. 

The plan of the paper is the following: in Section 2 we 
present the data and the definition of our samples, in Section 
3 we compare the redshift distribution of the radio galaxies sub- 
sample with the sample of all galaxies, in Section 4 we present 
the results of the color distributions and in Section 5 the bivariate 
luminosity functions. The radio-optical ratio distribution for our 
sample is presented in Section 6, while the stellar mass and star 
formation rate distributions are shown in Section 7. We estimate 
in Section 8 the star formation density evolution and in Section 
9 we present the conclusions. 

Throughout the paper we adopt a flat fl,„ - 0.3 and = 0.7 
cosmology, with H„ - 70 km s"' Mpc"'. Absolute magnitudes 
are given in the AB system and are expressed in the five stan- 
dard bands U (Bessel), B and V (Johnson), R and I (Cousins). 



Errors are estimated as poissonian fluctuation and refer to Icr 
confidence level. 



2. The data 

The radio data were obtained with the Very Large Array (VLA) 
in B configuration: a one square degree mosaic map with an ap- 
proximately uniform noise of ~ ILSfiiy (Icr) and with a 6 x 
6 arcsec FWHM gaussian resolution beam has been obtained. 
This map (centered at aaooo = 02:26:00 ^2000 = -04:30:00) has 
been used to extract a catalogue of 1054 radio sources complete 
to a limit of 80;uJy. A detailed description of the radio observa- 
tions, data reducti on, sources ext raction and radio source counts 
is reported in iBondi et al.l (l2003l) . while the radio-optical corre- 
lation (done with a likelihood ratio technique) and the proper- 
ties of radio objects with optical counterparts are presented in 
Cilie gi et"aLl ([2005 ). who found 718 radio sources with an opti- 
cal counterpart brighter than Iab - 25. The percentage of optical 
identifications of the radio sources varies from ~ 70%, for fluxes 
between 0.08 and 0.2 mJy, to ~ 65%, for fluxes between 0.2 and 
0.1 mJy, and decreases to 58% for brighter sources. Moreover, 
Bondi et al. (2007) performed a 610 MHz survey with the Giant 
Meterwave Radio Telescope (GMRT), finding that the median 
spectral index of faint radio sources (below 0.5 mJy at 1 .4 GHz) 
is significantly flatter than that of brighter sources. This fact is 
probably due to a relevant contribution below 0.5 mJy from a 
population of flat spectrum low luminosity compact AGN and 
radio quiet QSO. 

The optical data are part of both the F02 deep field of the 
VIMOS-VLT Deep Survey (VVDS) and the Dl field of the 
MEGACAM CFHT Legacy Surv ey (CFHT-LS). The firs t set 
comprises the B, V, R and I bands dMcCracken et al.Ll2003l) . ob- 
tained with the CFHT wide-field 12 K mosaic camera on the 
entire field. The completeness limits are ~ 26 for the B, V and R 
bands and ~ 25 for the I band. 

The second set c omprises the u*,g',r',i',z' bands 
(iMcCracken et al.L l2008 l) and covers an area of 1 deg^ 
centered on the VVDS F02 deep field; the completeness limits 
are 26.0, 25.8, 24.3, 25.1, 24.5, respectively. More over, we used 
the U KIDSS J and K data ( Umited at K~ 18, IWarren et all 
I2007h . covering the entire area. 

Starting from the VVDS photometric catalogue, a deep spec- 
troscopic survey has b een conducted for galaxies with Iab ^24 
(iLe Fevre et al.L l2"005h . with a sampling rate of ~ 33% over an 
area of 0.6 deg^ and with a spectral resolution of ~ 33A at 
7500A. Due to this sampling rate and the area covered, only 53 
radio sources (among the 718 with optical counterpart) have a 
spectroscopic redshift. Therefore, in order to increase the statis- 
tics, we decided to use the VVDS photometric redshift sample. 
The VVDS photometric redshift catalogue, the method to es- 
t imate redshifts and their quality are presented in Ilbe rt et all 
(2006). The typical error is o-(Az/(l + z)) ~ 0.029 and the re- 
liability is maximum in the redshift range [0.2 - 1.5]. We anal- 
ysed the differences between the photometric and spectroscopic 
redshifts for the 28 galaxies in the Radio Sample (see below), 
finding that the zero point is consistent with zero (at 0.05cr) and 
the error is 0.034. Therefore the properties of the photometric 
redshifts of radio galaxies are consistent with those of the opti- 
cal sample. 

Absolute magnitudes are computed following the method 
described in the Appendix of Ilbert et al. (2005): the K- 
correction is computed using a set of templates and all the pho- 
tometric information available. However, in order to reduce the 
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Table 1. Numbers of objects in our samples. The low/high redshift bins are [0-0.5]/[0.5-l.l] for all galaxies and type 3+4 galaxies 
samples and [0-0.7]/[0.7-1.0] for type 1+2 galaxies sample. 



redshift interval 




Complete Sample 




Control Sample 






all 


type 1 (z < 1.0) type 1+2 


type 3+4 


all type 1 (z < 1.0) type 1+2 


type 3+4 


0<z< 1.1 


430 


149 303 


127 


17880 3775 7414 


10466 


low z 


238 


176 


62 






high z 


150 


91 


59 







template dependency, the rest frame absolute magnitude in each 
band is derived using the apparent magnitude from the closest 
observed band, shifted at the redshift of the galaxy. With this 
method, the applied K-correction is as s mall as possible. 

Following the procedure adopted in IZucca etal] ( l2006h . for 
each galaxy the rest frame magnitudes were matched with 
an empirical set of SEDs, composed of four observed spec- 
tra (CWW iColeman et al., 1980 ) and two starb urst SEDs com- 
puted with GISSEL ( Bruzual & Charlotl 11993 ). The match is 
performed minimizing a variable on these templates at the 
redshift of each galaxy. Note that this method for the galaxy clas- 
sification is more refined than a simple color division because of 
the siultaneous use of all the available photometric data, which 
minimizes the influence of errors in a single band. 

Galaxies have been divided in four types, corresponding to 
the E/SO template (type 1), early spiral template (type 2), late 
spiral template (type 3) and irregul ar template (type 4); t ype 4 
includes also the starburst templates. iFranzetti et al.l(l2007l) have 
shown that there is a good correspondence between this classifi- 
cation and that obtained from spectral features, like e.g. the 4000 
A break up to z ~ lA. 

As shown by IZucca etal] (l2006h studying the luminosity 
function of these photometric types for the spectroscopic survey, 
type 1 (coiTesponding to early type galaxies) and type 4 galax- 
ies (corresponding to star forming late blue galaxies) show the 
most extreme behaviour regarding the evolution of the luminos- 
ity function, with the first class showing little evolution and the 
latter showing a strong evolution. The other two classes (2 and 
3) could be regarded as intermediate types in the sense that their 
properties are less extreme. 

The accurate knowledge of the behaviour of these classes led 
us to use the same classification: here, we assume that the results 
for the spectroscopic subsample holds also for the entire, photo- 
metric redshift based sample. This assumption is confirmed by 
the compaiison of the luminosity functions from the two datasets 
(O.Ilbert, private communication). 

From the 718 original radio-optical identifications (the 
"Original Sample"), we consider only objects with Mb < -20 
in the redshift interval [0-1.1] ([0-1.0] when considering type 1 
alone), where we have reliable photometric redshifts and well 
studied optical luminosity functions. We will refer to this sub- 
sample as the "Complete Sample". The magnitude limit assures 
that the considered galaxies are visible within the full adopted 
redshift range and therefore our sample has the same proper- 
ties as that of a volume limited sample. The different redshift 
limit for type 1 galaxies is due their higher K-correction, which 
influences their ab solute magnitude limit (see the discussion of 

IZucca et aL lHool)- 

The sample of aU galaxies from flie CFHT-LSA^VDS cata- 
logue has been cut at the same limits obtaining in this way the 
"Control Sample". In the upper panel of Figure [1] we show the 
apparent I-band magnitude distribution of the Original Sample 
(empty histogram) and of the Complete Sample (shaded his- 



togram). This histogram allows us to define a further cut in 
the Control and Complete Samples at Iab <24. This makes the 
two samples cleaner because at fainter magnitudes the estimated 
photometric redshifts are less precise and the fraction of wrong 
redshift increases. Moreover, this Umit in apparent magnitude is 
the same as for VVDS spectroscopic survey and therefore we 
are in the same framework as Zucca et al. (2006) work on the 
luminosity functions. 

Within these limits there are 430 radio sources. It is clear 
from the figure that the majority of the objects not included in the 
Complete Sample have magnitudes Iab >22. Among the radio 
sources not selected by our limits, 133 have Mb > -20 and 154 
match our absolute magnitude criterion but have z > 1.1, while 
the 35 Iab< 18 objects are underluminous to be included in the 
Complete Sample. 

Considering the types distribution within the Complete 
Sample, 169 objects are type 1 (149 within z < 1.0), 134 are 
type 2, while 74 and 53 are type 3 and type 4, respectively. Given 
the relatively low number of type 4 galaxies, we will consider as 
a single class (hereafter type 3+4) the late type galaxies. The 
Control Sample is formed by 4332 type 1 (3775 z < 1.0), 3082 
type 2 and 10466 type 3+4 objects (see Table [T]!. Therefore in 
our sample the percentage of radio emitting sources is ~ 4%, 
~ 3% and ~ 1% for the three considered classes, respectively. 
Here we assume that the emission mechanism of type 1 and type 
2 radio galaxies is mai nly due to AGN ( thus selecting the "ab- 
sorption line AGN" of [Best et all l2005l) . while for type 3 and 
type 4 it is induced by star formation. 

This assumption i s reasonab l e in t he light of the classifi- 
cation in Figure 9 of iBest et al.l (l2005h . which shows the plot 
logLradio/Mt versus Z)„(4000), where Lradio is the radio lumi- 
nosity, M, is the stellar mass and D„(4000) is the measure of the 
4000A break. In this plot there are two well separated clumps, 
not depending on logL,-adio/Mt, one at D„(4000) ~ 1.25 popu- 
lated almost completely by star forming galaxies, and the other 
at D„(4000) ~ 2.0 dominated by galaxies with AGN induced 
emission and LINERS. For the 31 radio galaxies in our sam- 
ple wdthjneasuredD,, (4000) our classification agrees with that 
of lBest etalJ (l2005l) . Considering the whole sample, our radio 
sources have almost all (~ 93%) a ratio 1 1.5 <logLradio/Mt < 14 
and this implies that the division between star forming and AGN 
induced radio emission spans the range [1.3 - 1.6] in D„(4000): 
this value roughly corresponds to the division between type 1+2 
and type 3+4 for the Control Sample. 

Another way to test our classification by looking at the 
upper panel of Figure [TT] (see Section 7): all radio galaxies 
with star formation rates (as estimated in the optical) less than 
1 - IM^yrs^^ (all of type 1+2) are far from the correlation be- 
tween SFR(optical) and SFR(radio) seen at higher rates, mean- 
ing that the AGN component of the radio luminosity is dominant. 
Moreover, we checked o ur classification w ith the spectra of the 
zCOSMOS lOK sample dLillv et all 120071) . by limiting the data 
set at Mb < -20 and z < 0.9, due to the different depth of the 
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survey. It results that only 6% of the type 1+2 objec ts falls un- 
der the dividing line (i.e. classified as star forming) of lBest et alj 
(2005) and 30 % of type 3+4 are just above the dividing line 
dBardelli et all 120091) . 

In the lower panel of Figure[T]we show the 1 .4 GHz flux his- 
togram of the Original Sample and the Complete Sample, while 
the long dashed line corresponds to type 1 galaxies and the dot- 
dashed line corresponds to type 3+4 objects. There is no major 
trend of types with flux: the fraction of type 1 galaxies decreases 
from 46% in the flux range [0.2-0.45] mJy to 32% in the range 
[0.08-0.2], while type 3+4 galaxies increases from 22% to 31%. 
In our case, no dramatic increasing percentage of star-forming 
galaxies is present at low fluxes, diffe r ently from early results 
(see the discussion in iGruppioni et all Il999l) . This conclusion 
would be similar if we adopt a fainter I magnitude limit or if we 
relax all our restrictions, in particular that at high redshift: in this 
latter case type 3+4 increases from 37% to 41%. Only the type 1 
show a decrease from 34% to 25%, possibly due to the blueing 
of ellipticals at higher redshifts, which could sh ift radio galax- 
ies fro m type 1 to type 2. However, as shown bv'Sevmo ur et alj 
(120081) . our adopted flux limit is just above the upturn of the star 
forming galaxies counts. 

2.1. Contribution from optically identified AGN 

In the optical spectroscopy, the emission of an active nucleus is 
revealed by the presence of typical features of emission lines, 
leading to the definition of two broad classes: the Broad Line 
AGN (BLAGN) and Narrow Line AGN (NLAGN). These AGN 
are known to be sometimes radio loud. These objects are ex- 
pected to be relatively rare within our limits and have spectral en- 
ergy distributions which could be misclassified within our clas- 
sification scheme (AGN induced radio emission associated to 
early type galaxies and star formation induced emission to late 
type galaxies). For this reason we consider as "contamination" 
the presence of these objects in both the Complete and Control 
samples. 

Here, we estimate how large the possible contamination in 
our sample by these objects is: unfortunately, a direct estimation 
is not possible because of the lack of spectroscopy for most of 
our radio sample. 

The presence of optical AG N in the spect r oscopi c 
VVDS survey h as been studied by Bongior no et al.l (l2007l) . 
iGavignaud et al.l 12006) and Paltani et al. (2008*), finding a sam- 
ple of 333 NLAGN and 63 BLAGN. Within our limits in redshift 
and absolute magnitude the ratio between type 1+2 objects ver- 
sus type 3+4 is 1/6 for BLAGN and 1/3 for NLAGN. Among 
these objects, 4 BLAGN and 2 NLAGN have radio counterparts 
in the Original Sample and 2 BLAGN and 2 NLAGN are part of 
our Complete Sample. 

The NLAGN optical counterparts have been assigned to the 
type 1+2 sample, while the BLAGN to the type 3+4 sample. 
The radio luminosities are logP(W Hz ') - 23.19 and logP(W 
Hz-i)= 22.52 for NLAGN and logP(W Hz-i)= 23.84 and 
logP(W Hz ')= 24.54 for BLAGN. These numbers come from 
the common area between the CFHT-LS survey and the spectro- 
scopic VVDS survey, which is approximately 28%. Moreover 
we corrected for the VVDS survey sampling rate (33%). Under 
the assumption that the contaminaton estimate obtained in this 
sub-area is representative of the entire area sampled, we expect 
in our sample a number of 21 ± 15 AGN objects in the power 
range logP(W Hz"') =23 - 24. 

As an independent way to estimate the contamina- 
tion, we correlated the positions of SWIRE survey sources 



dLonsdale et al.L l2b03h with our Complete Sample. For the sam- 
ple containing type 3+4 we have 85 common galaxies and a clear 
relation between redshift and the ratio of fluxes in the 24 //m 
and radio bands. Only 9 objects are significantly shifted towards 
lower values of the ratio, indicating the possible presence of an 
AGN. If the ratio between all galaxies and AGN is representa- 
tive, we expect a contamination in the sample of 15 ± 5 AGN, 
well consistent with the previous result. This approach is not 
possible for early type galaxies because the AGN and star form- 
ing populations are not separable. 

Considering the color distribution of all NLAGN, we find 
that for type 1 objects the B-I colors are more similar to the radio 
emitting objects, while the U-B colors have a flatter distribution, 
with a small peak at U-B~ 0.7. Also for type 4 galaxies, colors 
distributions have a similar behaviour. 

3. The Redshift Distributions 

In Figures |2] and [3] the redshift distribution of the radio galax- 
ies of our Complete Sample (solid lines) is compared with the 
Control Sample (dashed lines) for all galaxy types and for type 
1 and 3+4. Although not used in the following analysis, we show 
in the figures also objects with redshift larger than 1.1. Note that 
the Control Sample histograms have been rescaled by the ratio 
of the number of objects in the two samples. In principle, the 
two histograms are not expected to be exactly the same, being 
the radio sample cut both in optical magnitude and in radio flux 
while the Control Sample has only the first limit. The most im- 
portant peak in both the Control and in the Complete Sample is 
at z ~ 0.9 and is mostly due to early type galaxies. This max- 
imum is the result a the combination of t he presence of a huge 
large-scale structure at that redshift (see iLe Fevre et all l2005h 
and the increasing incompleteness of our Complete Sample be- 
yond z = 1 . 1 (or z = 1 .0 for type 1 galaxies) induced by our 
absolute magnitude limit. At lower redshifts, the histogram of 
the Control Sample is determined only by the increasing volume 
sampled. 

The other significant bump visible in the Complete Sample 
is centered at z ~ 0.3 and is mostly due to the radio flux limit. 
In fact, within our limits, the radio sample is complete for ra- 
dio powers of logP(W Hz ')>22 and for z < 0.2. At lower 
redshifts, the distribution is determined only by the increasing 
volume, while at higher redshifts the histogram is a combina- 
tion of the losses due to the radio flux limit and the increasing 
numbers of optical galaxies. However, also after a visual check, 
no obvious a - 6 segregation of radio galaxies is present corre- 
sponding to this bump. A clump at the same redshift is present 
in the Control Sample and in this case we detected the presence 
of a cluster (detected as a two-dimensional galaxy distribution) 
at aiooo = 2''24'"32^ and c^iooo = -4''14'54". Two radio sources 
have redshifts consistent with that of the cluster, one being at its 
centre. 

Finally, another difference between Control Samples and 
Complete Samples is an excess of radio galaxies at z ~ 0.6, 
which is the only one related to a structure, because of its nar- 
row redshift width. By plotting the a - 5 distribution of early 
type radio galaxies (which are expected to better delineate clus- 
ters and large scale structures) in the [0.5 - 0.65] redshift range, 
it is clear that they are mainly in the lower left (South-East) 
quadrant (Figure|4|l. By plotting the adaptively smoothed isoden- 
sity contours of all galaxies in this bin, we found that the over- 
density coincides with two galaxy clumps close to each other 
(~ 4 Mpc). This fact is consistent with the claim of Owen et al.l 
((19991) that cluster merging triggers radio emission (see also 
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Fig. 2. Redshift histogram of the radio galaxies in the Complete 
Sample (solid line) compared with galaxies in the Control 
Sample (dashed line). The Control Sample histogram has been 
rescaled by the ratio of the numbers of objects in the two sam- 
ples. 



radio source sample is of the same order because the correla- 
tion functions in the Control and Radio samples are similar (see 
Section 3.1). 

3.1. Correlation function 

An open problem concerning the properties of the radio galaxies 
regards their clustering properties. For this reason we computed 
the angular correlation function for all galaxies and for the radio 
loud objects divided in type 1 and type 3+4 galaxies. We did not 
detect any significant difference among the various populations; 
there is only a weak signal considering the bright and faint type 1 
radio galaxies, the ones fainter than logP(W Hz ')< 23.5 being 
more clustered than the brighter ones. 



Fig. 1. Upper panel: Shaded histogram: distribution of I mag- 
nitude for radio galaxies with Mb < -20 and Zpiwt < 1-1 (the 
Complete Sample). Empty histogram: all identified radio galax- 
ies (the Original Sample). Lower Panel: Dotted histogram cor- 
responds to the radio flux for all radio galaxies (the Original 
Sample), solid line represents the Complete Sample, dot-dashed 
line the type 3-1-4 objects, long dashed line the type 1 objects. 



iGiacintucci et al.L |2004|) . Unfortunately, no velocity dispersion 
estimate is possible because the two clumps lie just outside the 
region covered with the spectroscopic survey. 

Finally, note that the effect of cosmic variance on these plots 
is of the order of 15%, as estimated with similar samples, the 
same data (CFHT-LS) a nd the same classification method by 
iMcCracken et al.1 (l2008l) and m ore specifically for the VVDS 
survey by IVergani et al. (l2008l) . The expected variance on the 



4. Color Distributions 

The advantage of using redshifts is that one can use the rest- 
frame colors, by-passing all ambiguities led by the different type 
dependent K-corrections in adopting observed colors. Moreover, 
in order to disentangle real effects due to the presence of the 
radio source from the evolution of the optical host galaxy, it is 
necessary to make a relative comparison between radio galaxies 
of a given type with their "parent" population (in our case the 
Control Sample). 

The rest-frame B-I color distribution of early type galax- 
ies (type 1) which show radio emission is significantly different 
from that of the parent sample of early type galaxies (K-S prob- 
ability to be extracted from the same population of ~ 10"^), the 
first showing an evident peak at B-I ~ 1.4 (see Figure |5] upper 
panel), the second being broader, with approximately the same 
mode. 
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Fig. 4. Isodensity contours are the adaptive kernel smoothed a-6 
distibution of early type galaxies in the redshift range [0.5-0.65]. 
Points show the positions of the optical counterparts of the radio 
sources in the same redshift bin. 
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Fig. 3. Upper panel: Same as Figure |2]but for type 1 galaxies. 
Lower panel: Same as upper panel but for type 3+4 galaxies. 



The color distribution difference between radio emitting and 
parent sample early type galaxies is present only at higher red- 
shifts (0.5 < z < 1.0), while the distributions are statistically 
indistinguishable for z < 0.5. Moreover, the colors of early type 
powerful radio galaxies (logP(W Hz"')> 23.5) are not different 
from those of the Control Sample, while a significant difference 
(K-S probability of ~ 10"") is present between the color dis- 
tribution of faint sources (logP(W Hz ')< 23.5) and that of all 
early type galaxies. In all these cases the difference is due to a 
more peaked and narrower distribution of the radio sources with 
respect to the general population of early type galaxies. The red- 
shift and luminosity effect are related to each other in the sense 
that the higher redshift bin is populated mainly by the most pow- 



erful radio galaxies, while at lower redshifts the dominant popu- 
lation is the faint one. 

The difference in color distribution is present, although 
marginally significant, also in the U-B distribution (Figure |5] 
lower panel), where the significance of the difference of the ra- 
dio faint and bright galaxies with the Control Sample is 98% and 
95%, respectively. 

Considering the relative color distribution inside the radio 
sample, powerful radio sources tend to be bluer in B-I and redder 
in U-B than faint radio sources. This effect seems to be due to a 
relatively lower B band flux for radio faint objects. 

The difference in the B-I colors between radio loud and 
Control Sample is much more evident for type 3-1-4 objects, with 
a strong segregation of radio sources toward red colors with re- 
spect to the Control Sample (Figure |6] upper panel): this differ- 
ence is present, altough with far less significance, also in the 
U-B colors (KS probability of 0.01). In this case, the variable 
which drives this difference seems to be the redshift, as galaxies 
with redshifts below 0.5 have a less different distribution (KS 
probability of 0.01) than those at higher redshift (KS probabil- 
ity ~ 10 '^). This is due to the fact that at higher redshift the 
Control Sample containss more blue galaxies in the B-I colors 
than the Complete Sample. Opposite behaviour is detected for 
the U-B color, where no difference between radio loud and par- 
ent sample is found at higher redshifts, while a difference with a 
K-S probability of ~ 0.0007 is present at lower redshifts. 

An interesting question is whether the color difference of 
type 3-1-4 radio loud galaxies is due to some interaction between 
radio emission and optical spectrum or simply because redder 
galaxies tend to be more luminous in optical band (due to 
the well known color-magnitude relation for galaxies) and, 
given a radio-optical ratio, also in the radio wavelengths. To 
investigate this point, we took the Control Sample of type 3-1-4 
galaxies (10466 objects) and randomly assign to each galaxy 
a radio-optical ratio extracted from the observed distribution 
(see Section 6). With this value and given the optical luminosity 
of the chosen galaxy, we computed its radio power. Then, 



S.Bardelli et al.: The VVDS-VLA Deep Field: 



7 




U-B 

Fig. 5. Upper Panel: B-I color distribution of the radio emitting 
early type galaxies (solid line) compared with the same distribu- 
tion of the Control Sample (dashed line). The latter has been 
rescaled to take into account the different number of objects. 
Lower panel: U-B color distribution of the radio emitting early 
type galaxies and of the Control Sample. Symbols are the same 
as in the upper panel. 



we removed all objects whose simulated flux is below that 
of the VVDS-VLA radio survey and verified that the power 
distribution is similar to that observed. The resulting number of 
expected radio sources is 140 compared to the observed 143. 
Note that with this procedure we assumed the null hypothesis 
that the radio power is independent from the galaxy color but 
is depending only to the optical luminosity. As a result, the 
fraction of type 3h-4 mock galaxies with B-I< 0.5 is ~ 26%, 
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Fig. 6. Same as Figure|5]but for type 3h-4 objects. 

significantly higher than the observed ~ 8%. This means that, 
at a given optical luminosity, only the redder late type galaxies 
tend to have higher radio power, possibly indicating a larger 
amount of dust. In other words, higher star formation rates 
(which make a late type galaxy detectable in the radio bands) 
seem to occur in objects richer in dust (which makes the galaxies 
redder in the optical colors). 



5. The Radio Luminosity Function 

Luminosity functions are a powerful tool to investigate the evo- 
lution of astronomical objects. However, in the case of radio 
sources it is more correct to refer to a bivariate luminosity func- 
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tion, i.e. the composite distribution of optical and radio lumi- 
nosities. Therefore, one has to face the possible physical evolu- 
tion of the sources and biases induced by the flux limits in the 
two bands. This latter problem is particularl y severe in the opti - 
cal bands and at high redshift. As shown bv lllbert et al.l (|2004|) . 
in a magnitude limited sample a purely K-correction eff'ect can 
induce a variation of the population mix of red and blue galax- 
ies, even in the absence of evolution. Moreover, at different red- 
shifts any single photometric band records different regions of 
the galaxy spectra making the interpretation of the results more 
difficult. 

In order to overcome these difficulties, we choose to limit our 
optical sample to a B rest frame luminosity, i.e a quantity inde- 
pendent from redshift. The Mb - -20 limit permits a galaxy to 
be obseved in the whole adopted redshift range eliminating the 
optical constraint in the estimation of the observable maximum 
volume. The luminosit y function has been estimated with a stan- 
dard 1 / V„ax method ( S chmidtl 1 19681) in the redshift range [0- 1 . 1 ] 
([0-1.0] for early type galaxies). The radio K-correction has been 
estimated assuming that the radio spectrum is a power law (de- 
fined as 5 oc v") with a slope a - -0.7. As a check we recom- 
puted the luminosity function estimates with the observed 610 
MHz- 1.4 GHz indexes (for 310 objects within our limits) pre- 
sented by Bondi et al. (2007): the results remain unchanged. We 
present the results without consi dering the co rnpleteness correc- 
tion as determined in Table 3 of iBondi et alj (l2003i) . which sta- 
tistically corrects for radio sources lost at low fluxes. However, 
we have repeated the analysis also by weighting the l/V,„a.K of 
each radio galaxy by the corresponding flux-dependent correc- 
tion factor: the resulting corrections have the effect of an aver- 
age increase of the luminosity functions of ~ 10% for logP(W 
Hz"')< 23 in each redshift bin. For higher powers the correc- 
tions are larger (between 20 and 60%) but are in the same direc- 
tions for all redshift bins and therefore our conclusions do not 
change. 

In Figure |2] the radio luminosity functions of the Complete 
Sample (upper left panel), of early type (upper right panel ) and 
late type (lower right panel) galaxies are compared with the cor- 
responding low redshift results from iBest et af] (l2005l) (dashed 
lines). In particular, in the early type panel the red curve corre- 
ponds to type 1 galaxies, while the green curve represents the 
combined type 1 4- 2 sample, used to increase the statistics; in 
the late type panel the blue curve represents the combined type 
3+4 sample. 

The Complete Sample and the early type luminosity func- 
tions are in reasonable agreement with the lBest et al.l(l2005h esti- 
mate, also considering that their optical classification is slightly 
different, i.e. based on lines present in the galaxy spectra. The 
small systematic shift between our type 1 luminosity function 
and the corresponding iBest et af] (l2005h curve could be due to 
the different limiting optical magnitude, as the SDSS sample 
Umited to < -19. 

In order to check the effect of different optical luminosity 
limits, we derived luminosity functions for the Complete Sample 
with two different absolute magnitude limits, Mb > -19 and 
Mb > -2 1 . For this test we restricted the redshift range to z < 0.7 
in order to maintain the full visibility of the optical galaxies. We 
find that the two luminosity functions are essentially the same 
for logP(W Hz~')> 22.7, with the optically bright radio lumi- 
nosity function loosing ~ 20% of objects in the range of logP(W 
Hz ') [22-22.7]. Therefore we can conclude that all differences 
between our estimates and those o f iBest et all ( l2005b are due to 
the redshift evolution of radio sources. 



The type 3-1-4 galaxies luminosity function is flatter than the 
corresponding luminosity function of star forming galaxies of 
Best et al. (2005) and this remains true also considering only the 
most extreme population (type 4 galaxies). 

Assuming a power law shape for the radio luminosity func- 
tions, we find a slope of -1.78;^q q^ for type 1+2 galaxies and of 
-1.84+q [j for type 3+4 objects. 

In Figure |7] lower left panel, we show the comparison be- 
tween the optical luminosity function of radio emitting galaxies 
and the corresponding luminosity function of all galaxies in the 
same redshift bin (obtained following Zucca et al., 2006). Black 
curves correspond to all galaxies, red curves correspond to type 
1 objects, and blue to type 3+4 galaxies. Connected points are 
the 1/Vmax estimates for the radio samples, while solid curves 
are Schechter luminosity functions for optical galaxies, obtained 
with the standard paramet ric STY method for estimating op - 
tical luminosity functions dSandagei iTammann & Yahil 1979h . 
For clarity, the curves for various samples have been shifted ver- 
tically, by a value of +2 for type 1 and -1 for type 3+4. 

Working with luminosity functions means that the ratio be- 
tween the two curves for a given type in this panel is directly 
related to the probability for a galaxy to be radio emitter For lu- 
minous galaxies {Mb < -22.5) the probability is higher than for 
fainter ones and this seems to be independent from galaxy type. 
The probability for faint (Mb > -21) objects is ~ 7 x 10"- and 
1.7 X 10"^ for type 3+4 and type 1, respectively, and increases 
by more than a factor 3-4 for brighter galaxies. Note the flatten- 
ing of the rad io lumino sity functions at low powers: this effect 
is outlined by 'Cara&Xister (2008) in estimating the luminosity 
function with the l/Vmax in presence of evolving parent popula- 
tions. 

In Figure|8]we show the evolution of the luminosity function 
as a function of redshift. In the upper left panel the black curve 
represents radio galaxies of our Complete Sample in the [0.5- 
1.1] redshift interval, while the red one corresponds to the [0-0.5] 
bin. It is evident a significant evolution in the common range of 
radio luminosities reaching more than a factor 10 in the range 
logP(W Hz ') [24-25]. In the upper right and lower right panels 
the same is reported for type 1+2 and type 3+4 samples, respec- 
tively. We remind that given the different K-correction effect on 
absolute magnitudes, we limited the type 1 +2 sample to redshift 
< 1 .0. Moreover the limit between high and low redshift range 
is set to 0.7 for type 1+2 and to 0.5 for type 3+4 sample, in or- 
der to have reasonable numbers of galaxies per bin (see TablelTJ. 
These different limits have the advantage to show the maximum 
difference between samples. In order to increase the statistics, 
we show the type 1+2 sample instead of only type 1 galaxies. 
Considering the most extreme early type objects the results do 
not change but the significance decreases. The evolution of early 
type galaxies is mild, althought present, and it corresponds to a 
density increase of a factor ~ 5 in the logP(W Hz"^) range [23.5- 
24], while the evolution of type 3+4 sou rces is stro ng er. Thi s 
strong evolution is consistent wit h that of lAfon so et al. (2005J), 
which was found at z < 0.5, and ISmolcic et al., (2008) for the 
deeper VLA-COSMOS sample. Assuming an evolving law for 
luminosity and density evolution, of the form of (1 + zf and 
(1 + z)* respectively, we found that the best fit is yS = 2.76 
and k = 0.04 for typel+2 and /3 = 2.63 and k = 0.43 for 
type 3+4 objects. These values are consistent with the litera- 
ture dBrown et al 1 1200 It iHaarsma et al.L l2000t iHopkinsL l2004l) . 
However, the errors are quite large and the y6 and k parameters 
are strongly correlated: our aim is only to show that the early 
type galaxies present luminosity evolution with small density 
evolution, while late type galaxies show both evolutions. These 



S.Bardelli et al.: The VVDS-VLA Deep Field: 



9 



-3 



00 —4 
o 



-5 
-6 
-7 
-8 



1 Vj^ 1 1 1 1 1 1 


1 


fx 


All J 


1 1 1 1 1 1 1 1 


[ - 

\ : 

K , , , - 



22 24 26 28 
log(L,,,,„ [W Hz-i]) 




n 
O 

a 



CO 

G 
(U 
T3 

Xi 

B 
z 




ISO 
O 



CO 

e 

2 



22 24 26 28 
log(Lradio [W Hz-i]) 



■3 
■4 



-6 
-7 
-8 



1—^ — I — I — I — I — I — I — I — I — I — I — r 




J I I I I I 1_1J I I I I L 



22 24 26 28 
log(L,,,,o [W Hz-i]) 



Fig. 7. The radio luminosity functions of the Complete Sam ple for different galaxy types (points connected by solid lines) compared 
with NVSS/FIRST/SDSS survey results (iBest et al.Ll2005h . represented by a dashed line. Upper left panel: all galaxies. Upper right 
panel: early type galaxies; in red the type 1 sample and in green the type 1+2 sample. Lower right panel: late type galaxies; in blue: 
type 3+4 sample. Lower left panel: Optical luminosity function of radio emitting galaxies (connected points) compared with the 
luminosity function of the optical galaxies in the Control Sample in the same redshift bin (solid curves). Black: all galaxies; red: 
type 1 galaxies; blue: type 3+4 galaxies. The curves of the various samples have been vertically shifted for clarity (+2 for early type 
and -1 for type 1+2). 



results, obtained using a sample with a large redshift range and 
w ith high stat i stics, are c onsistent with the very early estimates 
of ICoUa et al] d 19751) and lAuriemma et all (Il977h . 

The detected evolution could have mainly two origins: a) 
the evolution of the corresponding optical luminosity functions, 
i.e. one has more objects and/or more optically luminous ones, 
which would translate in more luminous radio sources in the case 



of a constant radio-optical ratio; b) the evolution with redshift of 
the radio optical ratios; or a combination of both a) and b). 

In order to give an idea of the behaviour of the optical lu- 
minosity functions, which enter in point a), we show in Figure 
[8] (lower left pan el) the STY estimates (computed following 
iZucca et alll2006l) . for all galaxies (black), type 1+2 (green) and 
type 3+4 (blue) samples in the same redshift bins used for the 
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Fig. 8. Upper left panel: Radio luminosity function of all radio emitting galaxies with redshift below 0.5 (red curve) and in the range 
[0.5-1.1] (black curve). Upper right panel: the same but for type 1+2 sample and in the redshift intervals [0-0.7] and [0.7-1.0]. Blue 
dashed line is the de-evolved" luminosity functions (see text). Lower right panel: the same but for type 3-1-4 sample and in the 
redshift intervals [0-0.5] and [0.5-1.0]. Blue dashed line is the de-evolved" luminosity functions (see text). Lower left panel: optical 
luminosity functions for the Control Sample divided in the same types and in the same redshift intervals as previous panels. Dashed 
lines correspond to the high redshifts bins, while solid lines are the lower redshifts samples. Color codes are the following: black for 
the Complete Sample, green for type 1+2, and blue for type 3+4 samples. For clarity the curves of the variuos samples have been 
shifted in the magnitude axis {Mb + 1 for type 3+4 and Mb - 2 for type 1+2 samples). 



radio samples. The curves have been shifted for clarity in the 
magnitude axis by +1 for type 1+2 and by -2 for the type 3+4 
class. 

The main difference which can be seen in this panel is a 
brightening by one magnitude with a small steepening of the 
slope from low to high redhifts for type 3+4 objects. For type 



1+2 objects the main difference is that at high redshifts the 
galaxies are 0.2 mag brighter (from the best fit of M* in the 
Schechter function) and ~ 30% less numerous than today. These 
resul ts are consis t ent wi th those published for the VVDS survey 
(see IZuccaetaLL 120061 Table 3). 
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In the next section we will estimate the second parameter 
entering in the radio luminosity function evolution, i.e. the radio- 
optical ratio (which regards point b). 

We "de-evolved" the high redshift radio luminosity functions 
taking into account our conclusion regarding the optical lumi- 
nosity function and the radio-optical ratio for the two classes. 
This exercise aims to verify whether the detected evolution of the 
bivariate radio-optical luminosity function is the mirror of the 
evolution of both the optical luminosity function and the radio- 
optical ratio. 

In order to "de-evolve" the bivariate luminosity function of 
type 3 H-4 sample, we cut the higher redshift bin at Mg - -2 1 , as- 
suming an average one-magnitude evolution of the correspond- 
ing optical function (limited at Mb < -20 at z = 0) and no 
evolution of the radio-optical ratio. The luminosity function cut 
eliminates all optical galaxies entered at high redshifts in the 
radio sample and therefore responsible of the radio density evo- 
lution. 

For type 1-1-2 galaxies we corrected for density and luminos- 
ity evolution as detected for the optical luminosity function and 
a variation of ~ 1.4 in the radio-optical ratio (see Section 6). 

The resulting "de-evolved" luminosity functions are pre- 
sented as blue points in the right panels of Figure[8] The remain- 
ing discrepancy between these functions and the low redshift lu- 
minosity functions is likely due to contamination of BLAGN and 
NLAGN (see Section 2.1). 

Therefore, we estimated whether the difference between the 
"de-evolved" high redshift and the low redshift radio luminos- 
ity function of type 3-1-4 objects is consistent with the expected 
AGN contamination, assuming that all AGN are in the high red- 
shift bin. The number excess of the "de-evolved" luminosity 
function with respect to the low redshift one is 12, consistent 
again with the expected contamination of 15 ± 5 (see Section 
2.1). 

6. The Radio-Optical ratio 

In the previous section we showed the radio luminosity distri- 
bution for a sample with an optical absolute magnitude limit. A 
related way to show the relationship between the radio and op- 
tical emission is through the radio-optical ratio. Classically, it 
was estimated by using the quantity 5(,„yy) x iq("i-'^^-^)/^-^^ where 
S is the radio flux and m is the apparent magnitude of a radio 
galaxy. In our case, this formula is not satisfactory because it 
does not take into account the shift between the observed and 
the rest frame fluxes, which is important in our sample because 
of its large range in redshift. For this reason we compute the 
ratio between the radio luminosity {Lr) and the optical B band 
luminosity as 

Lr/Lb = U/lSm X 10"-4(5-48-M«) X 10^''] 

which represents a rest-frame radio-optical ratio, with both Lr 
and Lb in erg s However, this is not enough to estimate the 
correct radio-optical ratio distribution. In fact, the cut in radio 
power as a function of redshift induced by the flux limited sur- 
vey introduces a Malmquist bias in the ratio distribution, while 
this problem is not present for the optical data. For this reason, a 
given radio galaxy with an optical luminosity L'^ and radio lumi- 
nosity is visible to its maximum redshift z'„,g ^. determined by 
its radio luminosity and the limiting flux of the survey: the corre- 
sponding radio-optical ratio is therefore observable only within 
the corresponding V,',„^ . If one wants to compare the occurrence 
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Fig. 9. Frequency of the radio-optical Ratio in our sample. Solid 
line is the Complete Sample, dashed line corresponds to type 
1+2 galaxies and dotted line to type 3-1-4 galaxies 



of various radio-optical ratios within a fixed volume V, the radio- 
optical ratio of each single radio galaxy has to be weighted with 
the factor V/y,„a.v Moreover, comparing frequencies means that 
the histograms must be divided by the expected 2 W^ma v and 
not simply by the number of radio sources. 

The radio-optical distributions derived in this way for the 
Complete Sample and type 1+2 and type 3+4 samples are shown 
in Figure|9l The distributions are significantly broadened around 
the mean. For late type galaxies this is likely due to different dust 
content (seen as distribution in blue luminosity) corresponding 
to a given radio power, which is directly related to the star for- 
mation rate. For early type objects this is a measure of the AGN 
power variance, presumably related to the accretion efficiency, 
at a given optical luminosity. 

The radio-optical distribution of early type galaxies is shifted 
towards higher values with respect to that of late type objects. In 
order to study the evolution of the radio-optical ratio in the same 
redshift bins of the luminosity function, it is necessary to con- 
sider only those radio luminosity ranges which are always ob- 
servable in the considered redshift range. For this reason, we cut 
the type 1+2 sample at logP(W Hz ')> 23.1 and the type 3+4 
at logP(W Hz"')> 22.8, because the two classes have different 
high redshift limits. Although the number statistics are signifi- 
cantly reduced by these cuts, it is possible to say that no signifi- 
cant evolution is present for type 3+4 objects. This result is not 
surprising because if the shape of the initial mass function re- 
mains constant with redshift, the number of resulting radio emit- 
ting supernova remnants depends only on the total number of 
stars formed. However, we detected a slight dependence of the 
radio-optical ratio as a function of the optical magnitude, with 
the radio-optical distribution presenting a tail toward higher ra- 
tios for galaxies with Mb < -21.5 (FigurefTO] lower panel). We 
checked whether the increase of the radio-optical ratio for bright 
galaxies is due to the presence of an AGN using the upper panel 
of Figure[TT](see Section 7). The star formation values for bright 
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Fig. 10. Upper panel: frequency of radio-optical ratio as a func- 
tion of redshift for type 1-1-2 galaxies. Solid line represents low 
redshift (z < 0.7) galaxies, while dasehd line refers to high red- 
shift galaxies. To be comparable, the samples have been limited 
at logP(W Hz"')>23.1 Lower panel: the radio-optical ratio for 
type 3-h4 galaxies is plotted for objects with Mb < -21.5 (red 
dashed line) and Mb > -21.5 (blue solid line). 



galaxies are consistent with the logSFR(radio)- logSFR(optical) 
relation. 

For type 1 4-2 objects (Figure [TOl upper panel) the situation 
appears to be more complicated. It seems that there is a shift of 
a factor ~ 1.4 between the low redshift and high redshift bins. 
This would imply that the efficiency to feed the central black 
hole was a bit higher in the past. We used the reported factor 1 .4 



for deriving the "de-evolved" luminosity function described in 
Section 5. 



7. Star Formation Rate and Stellar Mass 

Star formation rates (SFR) and stellar masses for galaxies have 
been computed from the ph otometric quantitie s, following the 
method described in detail in lPozzetti et al.l (1200 7). by matching 
the data to a set of templates computed with various star for- 
mation histories and dust absorption. The star forma tion rates 
and th e stellar masses have been calculated using the 'Chabriej 
( 120031) initial mass function, while the equations used i n the lit- 
erature to compute SFR from the radio power use the ISalpetej 
(Il955h initial mass function. This implies that the SFR from the 
radio and optical band could be slightly different. For consis- 
tency, we converted our quantities by estimating the median dif- 
ference of\og(S¥Rsaipeier) " \og{S¥Rchahrier) ~ 0.186, when the 
two SFR are compared. In order to be consistent with other opti- 
cal papers, when we compare SFR and stellar mass distributions 
of radio loud galaxies and of the Control Sample we used the 
original Chabrieii (2003) initial mass function estimates. In the 
upper panel of Figure [TT] we show the star formation rate com- 
puted for the Complete Sample following this procedure based 
on the optical spectral distribution, compared wi th the star for- 
mation rate estimated following the formula of iHaarsma et al.l 
(l2000h (SFR(radio) = Li.4Gtf7/8.3610^°M o yr~' ) , which uses 
only radio powers at 1.4 GHz (see iCondom 1 1 992h . In the litera- 
ture there are other formulae to estimate the star fo rmation rate 
wit h different assumptions, as for example that of iBelll (|2003|) 
and ICram et a l. ( 1998) , which give rate s a factor 2.16 and 4.35 
lower than that of Haarsma et al.l ( l2000l) . respectively. 

The majority of early type galaxies (~ 70%) are segregated 
in the region with SFR(optical)< 2Mq this fraction in- 

creases to ~ 82% when considering only the extreme type 1 ob- 
jects. On the other hand, almost all (~ 86%) late type galaxies 
have star formation rates larger than 2Mq yr^\ This is not sur- 
prising because both the star formation rate and the type classifi- 
cation depend on the difference between red and blue-UV colors 
(mainly driven by the 4000 A break). The early type galaxies 
with SFR> 2Mq yr^^ show similar Mb with respect to the low 
star formation rate counterparts but a brighter ultraviolet color. 
This UV excess, which is responsible of the high star formation 
rate, has little, if any, impact onto the classification. 

There are some weak trends in the SFR(radio)-SFR(optical) 
plot for the two galaxy classes, type 1-1-2 showing an anticor- 
relation, while type 3-1-4 seem to be correlated. To evaluate the 
significance of the correlations we applied the Spearman rank 
test finding the values of pi+2 = -0.173 and P3+4 = 0.474 with 
a probability of the null hypothesis (no-correlation) of 0.002 in 
both cases. However, we assumed that the radio emission for 
type 1 +2 objects is AGN induced and therefore this relationship 
is to be regarded as radio power versus star formation rate. It is 
worth to note that the type 3-1-4 objects are not symmetrically 
distributed around the relation logSFR(radio)=logSFR(optical), 
being only ~ 19% of the points below the equality line. This 
result can be explained by the increasing role of dust with in- 
creasing star formation rates. Note that the factor between ra- 
dio power and star formation rate which maximize the frac- 
tion of objects w ith logSFR(radio)<logSFR(optical) is that of 
ICramet all d 19981) . leading to 34%. 

The lower pan el of Figure [TT] is analogous to Figure 3 of 
ICram et akl (Il998l) . The diagonal lines correspond to the charac- 
teristic time scales t= 10**, 10^ and 10'*' yr defined by Cram et al. 
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as MclSFRi Ar.H,. \Cmn\ et alj(ll998l) claimed that there is a ten- 
dency for galaxies that have already formed many stars to sup- 
port a higher current star formation; this effect is present also in 
our sample if we ignore the effect of different redshift. Dividing 
the sample in two redshift bins (above and below 0.5), we find 
that at the same stellar mass the star formation rate as estimated 
from the radio band is higher at higher redshifts. We checked this 
result also taking all radio galaxies with logP(W Hz"')> 22.8 in 
order to avoid biases due to different power ranges in different 
redshift ranges, finding the same result, i.e. there is no trend for 
the star formation rate with stellar mass, but the mean value of 
the star formation distribution was higher in the past. This means 
that the specific star formation rate is decrea sing with both red- 
shift and with stellar mass (see lVergani et al.L [2008). 

As done for the radio-optical ratio, in Figure [12] and [13] we 
plot the stellar masses and the star formation rate distribution 
weighted by V,„ax'- correcting stellar masses with the maximum 
volumes is equal to estimate the mass function. Our mass func- 
tions are not fitted by a Schechter function because the effect of 
increasing mass incompleteness with increasing redshift is im- 
portant in our large redshift ranges. However, we are interested 
in the comparison between radio and Control Sample and we 
assumed that the mass incompleteness is not a function of the 
radio properties and the same argument holds also for the star 
formation. 

Similarly to what has been found for the optical luminosity 
function of the radio emitting galaxies, also for stellar masses 
of type 1-1-2 objects the fraction of radio loud objects increases 
with the mass and almost all galaxies with M > 10" Mq are 
radio loud (upper panel of Figure [T2]l. 

In the lower panel of Figure[T2]we estimated the stellar mass 
functions in the two redshift bins we used for the luminosity 
functions. It is clear that the number density of radio sources 
is not changing with cosmic time, altho ugh the overa l l num - 
ber of massive galaxies is increasing (see [Pozzetti et all |2007|) . 
This probably means that the AGN is already active at high red- 
shift and no new activity formed in between. For what concerns 
the star formation rate, the distributions of radio loud and radio 
quiet type 1-1-2 objects are similar for SFR> 1 Mq yr"'. The two 
curves are only shifted by a factor less than 2. 

The distributions for late type galaxies are significantly dif- 
ferent, where the SFR of radio emitting objects is always lower 
than that of the Control Sample, although the difference is de- 
creasing with increasing star formation. This is true using the 
star formation rates estimated in the optical bands. When we use 
the star formation rate based on the radio power, we obtain the 
blue line of Figure [T3] more similar to the distribution of star 
formation rates of the Control Sample. Note that the behaviour 
shown for late type galaxies iri Figu re [T3] is similar to the rela- 
tions found bv lSevmour et al.l (120081) (see their figure 4). 

8. Star formation density evolution 

Assuming that radio emission from late type galaxies is entirely 
due to star formation, the estimation of the star formation den- 
sity evolution follows in a natural way. The most direct and used 
method is the determination of the radio luminosity function and 
its integral in different redshift bins. The integral of the radio 
luminosity function is directly related to the actual star forma- 
tion rate. The drawbacks of this method are essentially the poor 
statistics (few hundreds of radio objects at most are observed 
within the entire range of redshifts) and the fact that at different 
redshifts the sampled galaxies are different, being optically more 
luminous at higher redshifts. 
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Fig. 11. Upper Panel: Star formation rate computed from the 
optical versus the star formation rate determined from the ra- 
dio emission. Red and blue points correspond to type 1 -i-2 and 
3-1-4 objects, respectively. The vertical dashed line corresponds 
to the value of 2Mo3'ri while the solid line corresponds 
to logSFR(optical)=logSFR(radio) and is plotted for reference. 
Lower Panel: Stellar masses versus star formation rate derived 
from the radio emission for type 3-1-4 galaxies; squares represent 
the low redshift (z < 0.5) objects. 



Another way to approach the problem is that of considering 
the optical sample as a tracer of radio emission. In this case the 
star formation density could be written as 

SFRdensity^C f 0(L') ( f RD(R)L' dR] dL' 
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Fig. 12. Upper panel: Stellar mass functions of the type 1 +2 ra- 
dio sample (black solid line) and the type 1+2 control sample 
(red dashed curve). Lower panel: in red the low redshift ([0- 
0.7]) subsample for type 1+2 galaxies and in black the high red- 
shift sample ([0.7-1.0]). Dashed curves correspond to the control 
samples, solid curves to radio loud objects. 
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Fig. 13. Upper panel: Stellar mass functions of the type 3+4 ra- 
dio Sample (black solid line) and the type 3+4 Control Sample 
(red dashed line). Lower panel: Star formation rates distribu- 
tion of the radio rample (red solid line) and the Control Sample 
(black dashed line) for type 3+4 galaxies. The blue line corre- 
sponds to star formation rate computed from the radio emission. 



where R and D{R) are the radio-optical ratio and its probability 
distribution (derived as radio-optical distribution in Section 6), 
L is the blue absolute luminosity (see Section 2) , 0(L) is the 
optical luminosity function and C is the factor converting the 
radio lu minosity density to sta r formation density (assumed to 
be that of lHaarsma et all 120 00', see Section 7). 

Note that all the quantities in this equation have been esti- 
mated above. In particular, we take into account the difference 
between the optical luminosity function of the Control Sample 



and that of radio emitting galaxies and the dependance of the 
radio-optical ratio on the optical magnitude. 

We refitted the optical luminosity function in the B-band 
on the same sample and using the same method of IZucca et al.l 
(|2006|) but considering together type 3 and type 4 galaxies. For 
^{L) we used the fitted Schechter function, which depends on 
three parameters {a, M* and (p*). Errors are estimated consider- 
ing the variance in the fit of the optical luminosity function and 
considering the poissonian error on the radio-optical distibution 
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Fig. 14. Star formation rate density as a function of the red- 
ds hift. The red po i nts ar e our estimates, t he black points those 
of iHaarsma et all (|2000|) as c orrected by iHopkinsI 1200 4), the 
green solid hne is the fit of iHopkinsI (I2004|), the black dot- 
dashed li ne is from Cram et al. ( 1998). The blue squares are the 
points of ISmolcic et al., (2008) and crosses are the estimates of 
ISevmour et al. 120081) . Large black squares with error bars cor- 
respond to the VVDS dust uncorrected values, while red and 
blue dashed lines correspond to the star formation rate density 
estimated from near-IR and Ho- data as reported bv lTresse et alj 
(I2007h . 



function. In fact, we computed various Schechter functions on 
the border of the one sigma error ellipses and considered the 
maximum and minimum functions as part of the error induced 
by the optical luminosity function. Similarly, we have taken as 
upper and lower limit for the radio-optical distribution function 
the one sigma fluctuations in the histogram of Figure|9]and again 
we assumed the maximum and minimum of the computed star 
formation density as error induced by the radio-optical distribu- 
tion. Finally, we added in quadrature these values with the pre- 
vious ones. 

For the last two redshift bins, where the slope a of the optical 
luminosity function is not constrained, we fixed the value of the 
slope to Q- = -1 .33, which is the fitted value of the two bins con- 
sidered together. Therefore, for the last two bins the errors given 
by our procedure are underestimated. The values are reported 
in Table |2] In Figure [14] the estimated star formation density is 
compared with values found in the literature both from radio and 
other bands. The radi o data has been converted to the formula of 
iHaarsma et al. I (120001) . 

There is an overall consistency between all the estimators, 
even if our star formation density in the lowest redshift bin is 
higher than the other estimates. However, considering that the 
real error bars could become a bit larger by taking into account 
other sources of statistical variation than error in the optical lu- 
minosity function fit alone, our point at z = 0. 1 could be regarded 
as still consistent with othe r estimators. Moreover, also the first 
bin of ISmolcic et alj (l2008h . whose data in term of limiting ra- 



dio flux are similar to ours, is high. By comparin g our values 
with the VVDS dust uncorrected values found by Tresse et akl 
(2007) using ultraviolet data, we find that at all redshifts there is 
an approximately constant shift of A logSFR^^,.,,,,, = 0.73 + 0.02. 
There is also a weak indication of a ~ 10% increase of this cor- 
rection with redshift. N ote that the value is v ery similar to the 
correction of ~ 0.74 of Ho pkins et alj ( 1200 Ih (see equation 7) 
required for the FUV 1500A (for which the lTresse et"an.l2007l 
SFR densities are computed) and a star formation rate of logSFR 
~ 1. 

Moreover, our star formation density is fitted by an evolution 
of ( 1 + z)^ '^, consi s tent w ith the ultraviolet light density evolu- 
tion of iTresse et al] (l2007h . 

What is important here is the conclusion that the optical light 
could be used as tracer of the radio emission at high redshifts, 
provided that no significant variation of the radio-optical ratio 
is present. Note that we checked this point because our radio 
sample cover, even if with small statistics, a large fraction of the 
range over which the star formation history has been estimated. 

9. Conclusions 

In this paper we analyzed the rest frame properties of radio loud 
galaxies of the VVDS-VLA Deep Field survey. In order to avoid 
all difficulties intrinsic to samples containing both optical and 
radio limits, we have chosen to a) work using rest-frame quan- 
tities; b) limit the sample to an absolute magnitude Mb < -20 
and c) consider only objects with z < 1.1 (< 1.0 in the case of 
early type galaxies). The two last points correspond to having 
an optical volume limited sample, i.e. no loss of optical galaxies 
with redshift. Therefore, the only limit which remains is the ra- 
dio flux limit (~ SOfiJy). In this way we can control the optical 
properties of radio galaxies, although with paying the price of 
smaller statistic s and lo osing all optically faint objects. 

Following iZucca et al.l (120061) . we divided the galaxies in 
four spectrophotometric classes corresponding to early type 
(types 1 and 2) and late type (type 3 and 4) galaxies. 
Furthermore, we assumed that early type galaxies have AGN- 
induced radio emission and late type galaxies have star- 
formation induced emission. The aim of this work was to inves- 
tigate the possible evolution of radio galaxies knowing a priori 
the optical behaviour of all galaxies with our cuts in the entire 
CFHT-LS/VVDS survey (taken as Control Sample). 

From the redshift distribution of radio sources we detected an 
overdensity corresponding to two interacting clusters at redshift 
z ~ 0.6, while the other redshift peaks appear to be associated 
with a large scale structure (superclusters). It has been impossi- 
ble to investigate further this point because of the relatively small 
number of spectroscopic redshifts and the smoothing effect due 
to the use of photometric redshifts. The angular two points corre- 
lation function of radio loud galaxies does not show significant 
differences with respect to that of optical galaxies, suggesting 
that these objects are not living in the most extreme environ- 
ments. By estimating the distribution of blue luminosities, stel- 
lar masses and star formation rates, it results that the probability 
for a galaxy to be a radio emitter increases significantly at high 
values of these parameters. 

Early type (type 1) radio loud galaxies show some differ- 
ences to respect to their radio quiet counterparts in the colors 
distribution. In particular, radio faint galaxies show a peak in 
the B-I color at B-I~ 1.45, while radio bright galaxies show no 
difference with the Control Sample. This class of radio sources 
shows evolution in the bivariate radio-optical luminosity func- 
tion, mainly due to luminosity evolution. Studying the radio- 
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Table 2. Parameters of the optical luminosity function and value of the star formation rate. The symbol * means that errors are 
underestimated. 
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optical ratio distribution we show that for redshifts above 0.7 
at fixed blue luminosity the galaxies were more radio luminous. 
This effect is somehow compensated in the bivariate luminosity 
function because of the simultaneous decrease in the type 1 +2 
optical luminosity function. This luminosity evolution is likey 
due to a different efficiency to feed the central black hole and not 
to an increased contribution by star formation in the radio band. 
This could be seen also by considering the results of Figure [TT] 
(upper panel) divided in redshift bins: the fraction of galaxies 
with star formation rates (as computed in the optical band) lower 
than 2Mq yr"' is the same at high and low redshifts. This means 
that, unless dust absorption and star formation exactly compen- 
sate each other, no significantly increased star formation for this 
galaxy type is present. 

Type 3 and type 4 radio loud objects are significantly redder 
with respect to the control population. These galaxies show a 
significant evolution in the bivariate luminosity function both in 
luminosity and density: this behaviour can be explained if we 
take into account the one-magnitude evolution of this class of 
objects detected in the optical, without invoking large variations 
in the radio emission properties (like the radio-optical ratio). In 
fact, the only dependence of the radio-optical ratio is, in this 
case, based on the optical luminosity. Using the estimated radio- 
optical ratio we conclude that the difference in colors between 
radio loud and radio quiet objects is physical, and not induced 
by statistical effects: we speculated that radio loud star forming 
galaxies contain on average more dust. 

The star formation rate as computed from radio power and 
from optical colors shows a significant correlation si milar to the 
relat ion betwen radio and UV-based SFR shown in ICram et akl 
(1199 8). By plotting the star formation rate computed from radio 
luminosities and the stellar mass, we show that at higher red- 
shifts the population has higher instantaneous star formation at 
a fixed stellar mass. In other words, the specific star formation is 
decreasing with redshift. 

With the knowledge of the parameters for the radio loud type 
3 and type 4 galaxies, we show that it is possible to use opti- 
cal galaxies as tracers of the radio emission at various redshifts. 
This permits to derive an estimate of the star formation history, 
which resulted to be consistent with the results of other bands 
and methods. 

In summary, from our data, the AGN induced radio emission 
is increasing with redshift because the host galaxies increased 
their radio-optical ratios, with the AGN already in place at high 
redshifts. On the contrary, the evolution of late type radio sources 
is a direct consequence of the strong evolution seen the in optical 
band and is due to the evolution of the star formation rate with 
redshift. 
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